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Charmonium dissociation temperatures are studied in a quenched anisotropic lattice QCD with 
standard plaquette gauge action and 0{a) improved Wilson fermion action. Simulations are car- 
ried out at temperatures in the range 0.88 to 2.3 T^. From the meson correlators, we first 
subtract the contribution of constant mode, which was reported to mislead the analysis of, in par- 
ticular, P wave signals. We then calculate effective masses and Bethe-Salpeter wave functions 
for ground (IS, IP) and excited states (2S, 2P) using the multi-state variational analysis. To dis- 
tinguish between bound states and scattering states, we apply two methods: First, we compare 
effective masses for charmonium correlation functions with finite spatial momenta under differ- 
ent spatial boundary conditions. Since the scattering state energies are sensitive to the boundary 
conditions, we expect finite volume effects when the charmonium dissociates. Second, we study 
the Bethe-Salpeter wave function, which should become broad when the charmonium bound state 
turns into a scattering state. With both methods, we have fond no clear evidences of dissociation 
for the ground and exited charmonium states up to 2.3 Tc so far 
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1. Introduction 

Studying heavy quarkonia properties above Tc is important to understand how quark-gluon 
plasma (QGP) is formed in heavy ion collisions. In particular, dissociation of charmonia in de- 
confined phase is an interesting subject because //y/ suppression [|l|] is one of the most important 
signals of QGP. Recently, the sequential //y/ suppression scenario [^, in which dissociation of 
heavier charmonia {Xc, w') pl^Y important roles for the suppression of total 7/y/ yield, has been 
proposed to understand the experimental results in heavy ion collisions at SPS and RHIC. 

Lattice QCD studies of charmanium spectral function using the Maximum Entropy Method 
(MEM) suggested that S wave charmonia, such as r\c and J /^r, survive up to 1 .5 |, ^ |, ^ |8|] , 
while P wave charmonia, such as Xc^, dissociate just above [^, ^, ||]. 

Although these results look supporting the sequential 7/ suppression scenario, there are still 
some problems. First, the meson correlators receive sizable contribution of the constant mode 
above Tc, especially in P wave correlators!^. From an analysis in which the constant mode con- 
tribution was properly taken into account, the thermal effects in charmonium correlators including 
those of P waves turned out to be quite small up to 1.4 Tc. This is in contrast to the conclusion 
of previous MEM studies on P wave charmonia. Furthermore, choice of the default model for 
the spectral function is known to affect the MEM results rather sensitively. A crosscheck without 
Bayesian-type analyses is needed. 

Let us first discuss what we expect on the lattice when a charmonium dissociates. It should be 
noted that, on a lattice with finite extent, spectral functions consist of discrete spectra only, also at 
T > Tc. (These discrete spectra appear as broad peaks in approximate calculations such as MEM.) 
Below Tc, we expect to have discrete peaks corresponding to ground and exited bound states. When 
a charmonium bound state fully dissolves above Tc, the peak corresponding to the bound state will 
vanish and other peaks corresponding c-c scattering states may appear near the vanished peak. 
Furthermore, we expect that the wave function for the peak will have a localized shape when the 
peak corresponds to a bound state, while it will have a broad shape extending to large distances 
when the peak corresponds to a scattering state. 

In this paper, we report on our study of charmonium dissociation. In principle, when a high 
precision data of the correlation function is available up to large distances, the constant mode can be 
identified by conventional analyses too. In practice, however, with current accuracy and the range 
of sensible data, we think that it is safer to explicitly subtract out the contribution of the constant 
mode from the conelators. We apply the midpoint subtraction method developed in [^]. Avoiding 
Bayesian-type analyses, we then adopt the multi-state variational method to extract ground and 
exited state masses both in S and P waves [p^]. At the same time, we calculate the Bethe-Salpeter 
wave functions of these states from the spatial correlation function between c-c quarks [||]. 

To distinguish bound states of c-c quarks from their scattering states for the extracted states. 



we study the spatial boundary condition dependence of the energy spectrum at finite volume[|llp: 
The energy of scattering states depends on its relative momentum which is quantized according to 
the spatial size and boundary conditions. On the other hand the spectrum of the bound states does 
not change against such exchange of boundary conditions. The Bethe-Salpeter wave functions for 
a bound state will be compact, while those for scattering state will extend and will change its shape 
depending on the spatial lattice size. Combining these tests, we may examine if the charmonium is 
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Table 1: The parameter A,- of the smearing functions. 
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2. Multi-state variational analysis 

The charmonia correlation matrix in Euclidian space-time is defined by 

Co(O^I(of(^,Ooy(o,o)t). 



(2.1) 



with r = 75, t;-, 1 (/ = 1, 2, 3) for pseudo-scalar (Ps), vector (Ve) and scalar (Sc) channels, respec- 
tively. Here {Of | / = 1, 2, • • ■A'^state} are smeared meson operators defined by 



of =Y^(Oi{y)(Oi{z)q{x + y,t)Tq{x + z,t), 



(2.2) 



with Gaussian smearing functions (Oi{x) = exp(— A, |.xp). Table [T] shows our choice of the smearing 
parameters A,. 

According to [^, constant mode effects are large in the deconfined phase, because the constant 
mode is due to wraparound contributions of single quark propagators, which are suppressed in the 
confined phase. To separate out the constant mode contribution from meson correlators, we study 
midpoint subtracted correlators C{t) defined by 



C(0 = C{t)-C{Nt/2) 
= [cQ + Y^Ck cosh 

V k=i 



Nr 

m\t- — 



*:=0 



l^Ck sinh^ 



k=i 



2 I 2 



(2.3) 



Here C{t) is the matrix of C,y(f) defined in (g^, cq is a contribution of the constant mode and Nt 
is the temporal lattice size. 

We extract ground and excited states in C by solving the eigenvalue problem. 



C{t)vu = h{t;to)C{to)vk {k = l,2, A^state). 
Effective masses of the diagonalized states are given by 



(2.4) 
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(2.5) 
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Table 2: The number of gauge configurations. 
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Next, we define wave functions in terms of the Bethe-Salpeter (BS) amplitude: 

BSi{r,t) = £(^(x,0rs„k^(x + r,0Of-(0,0)^) 

^(i/Ai(r));cosh ' ' ^' 



(2.6) 



Here {^fk{r))i is the wave function for the k-\h state created by the smeared source operator 

To study S waves, we adopt Fgrc = ^&ak = 75 and yt for Ps and Ve charmonia, respectively. For Sc 

charmonium in P wave, we use a derivative source Fsrc = 7(5; — di) to enhance the signal, with 

Tsnk = Yi- 

To seprate out the constant mode, we define midpoint subtracted BS amplitudes by 

BS{r,t) =BS{r,t)-BS(f,Nt/2). (2.7) 
Finally, the diagonahzed wave functions are given by 

Y,iBSi{rQ,t)Vik 

where Vik = {vk)i is the i-th component of the eigen vector v^, and ro is the normalization point. 

3. Numerical results 
3.1 Lattice setup 

We study quenched QCD with the anisotropy ^ =as/at = 4. The gauge action is the standard 
plaquette action with j8 = 6.10 and a bare anisotropy parameter Jg = 3.2108. For quarks, we adopt 
0{a) improved Wilson fermion action with tadpole-improved tree-level clover coefficients. Our 



actions and their parameters are the same as those in [ |12| ] except for the choice of the Wilson 
parameter, r = 1. The spatial lattice spacing is a^^ = 2.030(13) GeV. Simulations are performed 
on *Ni lattice, with A^, = 16, 20 and 32 and Ni = 32, 26, 20, 16 and 12, which correspond 
to r = {NiUt)^^ = 0.88 Tc, 1.1 Tc, 1.4 Tc, 1.8 Tc and 2.3 Tc respectively. (Critical temperature 
corresponds to Nt ^ 28.) The number of gauge configurations is shown in Table Q We choose the 
Coulomb gauge fixing to calculate the wave functions. 

3.2 Effective masses 

In Figure [l] we compare the ground state's and the first excited state's effective masses of Ve 
and Sc charmonia as functions of T between different boundary conditions, i.e. periodic (PBC), 
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Figure 1: The effective masses for Ve (left) and Sc (right) charmonia. Lower and upper symbols correspond 
to the ground and first excited states, respectively. Red, black and blue symbols indicate the effective masses 
with PBC, APBC and MBC, respectively (see text for the definition of MBC). Short bars shown near the left 
end of each figure indicate the boundary condition dependence of energy spectra for the case of free quarks. 



anti-periodic (APBC) and a mixed (MBC) boundary conditions. Here MBC is defined by APBC in 
the x-direction combined with PBC in y and z-directions. At zero temperature, the ground state and 
the first excited state of Ve channel correspond to J/y and y' respectively, and the ground state of 
Sc channel corresponds to XcO- The effective masses are calculated by the multi-state variational 
analysis with 4x4 correlation matrix on the 20^ x Nt lattice. Here we choose smeared operators 
with the parameter A3 , A4 , A5 and A^ of Table |I] to obtain the largest overlap with the ground state 
and the first excited state. 

To see a typical magnitude of the boundary-condition dependence in the effective masses, we 
have studied the energy spectra for the case of free quark in a box of (2 fm)-', assuming that the 
charm quark mass is 1.3 GeV. Results are shown by short bars around the left end of each plot in 
Fig.|[ (Only relative differences in the vertical coordinates are relevant for the bars.) For ground 
state P waves, we obtain the largest gap between PBC and MBC. Therefore, we may expect a mass 
shift of about 200 MeV (500 MeV) for a ground (first exited) state when the charmonium fully 
dissolve. 

In Fig.|l], we find no evidence of mass shift with such magnitudes, for neither the ground states 
nor the first excited states, up to a quite high temperature of about 2.3 T^. The results for Ps 
charmonia are similar. This suggests that these charmonia bound state persist up to about 2.3 Tc. 

3.3 Wave functions 

To check the bound-state characteristics of the charmonium states, we study the wave functions 
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Figure 2: The wave functions for Ve (top) and Sc (bottom) charmonia at 2.3 Tq- Horizontal axis is the 
spatial distance in the it-direction. Left panels show the results for ground states and right panels show those 
of the first excited states. Red, green and blue symbols indicate data obtained on Ns = 32, 20 and 16 lattices. 
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corresponding to the mass eigenstates diagonalized by the multi-state variational analysis discussed 
in the previous subsection. When a charmonium dissociates, corresponding wave function should 
extend to large distances and show a sensitive dependence on the spatial size of the box. Figure 
^ shows the ground state's and the first excited state's wave functions defined by (2^) for Ve and 



Sc charmonia at 2.3 Tc, obtained on three different spatial lattices of 32^, 20^ and 16^. We find 
that, for all charmonia we study, the spatial size dependence is small, and the wave functions 
are compact. We thus find no sign of c-c scattering states up to 2.3 Tc, in accordance with the 
observation discussed in the previous subsection. 

4. Conclusions 

We investigated if charmonia dissociates in high temperature QCD on quenched anisotropic 
lattices. Adopting the multi-state variational analysis to extract ground and first excited states, we 
examined spatial boundary condition dependence of effective masses and investigated the shape 
and spatial size dependence of wave functions for these states. From these studies, we found no 
sign of scattering states nor clear evidence of dissociation, both for the ground and first exited states 
of charmonia in Ps, Ve and Sc channels up to a quite high temperature of 2.3 Tc, so far. 

We thank the members of the WHOT-QCD Collaboration for innumerable discussions and 
suggestions. This work is in part supported by Grants-in-Aid of the Japanese Ministry of Education, 
Culture, Sports, Science and Technology (Nos. 17340066 and 19549001). Numerical calculations 
were performed on supercomputers at RCNP, Osaka University. 
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